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We study the role played by multiply-connectedness in the time evolution of the energy E(t) of a 
radiating system that lies in static flat space-time manifolds M4 whose t = const spacelike sections 
M3 are compact in at least one spatial direction. The radiation reaction equation of the radiating 
source is derived for the case where M3 has any non-trivial flat topology, and an exact solution 
is obtained. We show that the behavior of the radiating energy E(t) changes remarkably from 
exponential damping, when the system lies in 1Z 3 , to a reverberation pattern (with discontinuities 
in the derivative E(t) and a set of relative minima and maxima) followed by a growth of E(t), when 
Mz is endowed with any one of the 17 multiply-connected flat topologies. It emerges from this 
result that the compactness in at least one spatial direction of Minkowski space-time is sufficient 
to induce this type of topological reverberation, making clear that topological fragilities can arise 
not only in the usual cosmological modelling, but also in ordinary flat space-time manifolds. An 
explicit solution of the radiation reaction equation for the case where AI3 = 1Z 2 x S 1 is discussed in 
details, and graphs which reveal how the energy varies with the time are presented and analyzed. 



I. INTRODUCTION 

In general relativity, as well as in any metrical theory of gravitation of some generality and scope, a common approach 
to cosmological modelling commences with a four-dimensional space-time manifold (representing the physical events) 
endowed with a Lorentzian metric (necessary to ensure the local validity of the well established special relativity 
theory) . Representing the physical phenomena in this manifold we have fields satisfying appropriate local differential 
equations (the physical laws). Finally, fields and geometry are coupled according to the corresponding gravitational 
theory one is dealing with. The space-time geometries arising as solutions of the gravitational field equations constrain 
to some extent the dynamical behavior of the physical fields. This eminently metrical approach to model the physical 
world has led a number of physicists to implicitly (or explicitly) restrict their studies to purely geometric features of 
space-time, either by ignoring the role of topology or by considering just a limited set of topological alternatives for 
the space-time manifold. 

If, on the one hand, the topological properties of a manifold antecede and are more fundamental than the metrical 
features and the differentiable structure on which tensor analysis is based upon, on the other hand, it is well known 
that geometry constrains but it does not dictate the topology of a space-time. It is therefore important to determine 
whether and (or) to what extent physical results concerning a space-time geometry depend on or are somehow 
constrained by the topology of the underlying manifold. 

In a recent work we studied the role played by topology in the time evolution of the energy E(t) of a radiating 
system in static flat FRW space-time manifolds A^4 whose t = const spacelike sections M3 are endowed with different 
topologies, namely the simply-connected Euclidean manifold and six topologically non-equivalent flat multiply- 
connected orientable compact 3-manifolds - fl. Clearly these flat space-time manifolds M.± are orientable and 
decomposable into A^4 = 7Z x M3 . The radiating system we have examined, which has been used to study a wide 
class of oscillating phenomena Q - flqj , is formed by a pointlike harmonic oscillator (energy source) coupled with 
a massless scalar field (scalar radiation waves propagating at speed of light). Through a numerical integration of 
the time evolution equation for the harmonic oscillator (with the topological constraints suitably considered), we 
have shown that there is an exponential damping of the energy E(t) of the harmonic oscillator when the spacelike 
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section M.3 is the Euclidean space 1Z 3 , whereas when M.3 is endowed with one of the six compact and orientable flat 
topologies the energy E(t) exhibits a loosely speaking reverberation pattern, with discontinuities in the derivative 
of E(t) and relative minima and maxima (both the discontinuities and the extrema are due to the travelling waves 
which "reflect" from the "topological walls"), followed by a growth of the energy with the time. We have also shown 
that, for these six compact cases, the energy E(t) diverges exponentially when t — > 00 , in striking contrast with the 
damping of E(t) in the 1Z 3 case. This unexpected divergent behavior of the energy for flat space-times with compact 
spacelike sections illustrates that (totally) compact topologies may give rise to rather important dynamic changes in 
the behavior of a physical system. This type of sensitivity has been referred to as topological fragility and can occur 
without violation of any local physical law fllTfl . 

In the present work we extend these investigations by performing a rigorous non-numerical study of the behavior 
of the same physical system under a less restrictive topological setting, namely we shall consider static flat space-time 
manifolds M.4 whose t = const spacelike sections M.3 are compact in at least one spatial direction. Moreover, no 
assumption will be made as to whether or not these 3-manifolds are orientable. This amounts to saying that the 
results of the present work hold for all flat space-time manifolds whose t — const spacelike sections are endowed with 
one of the 17 flat topologies discussed by Wolf Ellis ||, and others Q (see also Gomero pi). 

It should be stressed from the outset that in spite of being the same physical system discussed in Q , the present 
work generalizes the results of that article in several respects. Firstly, the underlying topological setting is now much 
more general in that here we do not assume that the spacelike 3-manifolds M3 are compact and orientable. Secondly, 
contrarily to the numerical integration performed in |lj , in the present paper we have obtained a closed exact solution 
of the evolution equation for the harmonic oscillator in the above-mentioned general topological setting. Thirdly, the 
case study we have examined in section [V, where we obtain the explicit solution for the radiation reaction equation 
and discuss in details the graphs which reveal the time evolutions of E(t) and E(t), is completely new. Finally here 
we have also performed a lengthy discussion on the conservation and on the balance of the total energy of our physical 
system, making apparent how (here as well as in ref. jp) the total energy of the whole system is conserved for any 
time t > 0. 

The plan of this article is as follows. We derive (section |l|) the time evolution equation for the harmonic oscillator 
(also referred to as radiation reaction equation) when M3 is endowed with a non-trivial flat topology, i.e., when the 
t = const spacelike sections are any multiply-connected (compact in at least one direction) flat manifold M3. It turns 
out that the evolution equation for this case is formally the same as that obtained by Bernui et al. 13] for the case 
where M3 is any (totally) compact and orientable 3-manifold. However, they differ in the non-homogeneous term 
that accounts for the difference in the topology — distinct degrees of connectedness of the spacelike sections give 
rise to different non-homogeneous terms. We also emphasize in section || that since our system is composed only of 
test fields in static flat space-time backgrounds, there is no need to consider the back-reaction of the fields on the 
geometry. 

We present in section III an exact solution of the evolution equation for the harmonic oscillator in flat space- 
time whose t = const spacelike sections are any multiply-connected flat 3-manifold. Through a heuristic analytical 
reasoning we also show in section [II that this exact solution Q(t) and the corresponding energy of the harmonic 
oscillator E(t) both exhibit a divergent asymptotical exponential behavior. Clearly the energy can be limitlessly 
"extracted" from the interaction term in the multiply-connected cases only because the system permits this extraction, 
since it is not bounded from below. But as we shall discuss in this work it is indeed the topology that "excites" this 
available "physical mode" of our system. The balance and conservation of the total energy of our system is also 
discussed in section III, where we show that it is finite and conserved for all finite time t > 0. 

Explicit exact solution of the evolution equation for the case where M3 = 7Z 2 x S 1 is discussed in details in 
section |v[ For this special case we also present and analyze graphs which reveal how the amplitude Q{t) and the 
energy E{t) of the overdamped harmonic oscillator vary with the time. The graph of E(t) presents a reverberation 
pattern with relative minima and maxima, and discontinuities of the derivative E(t) . The topological origin of the 
reverberations for the overdamped oscillator is also made clear and stressed in that section. The graphs in this section 
also confirm (within the limit of accuracy of the plots, of course) the divergent asymptotical exponential behavior 
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for both Q(t) and E(t), which is shown to take place in section III. These results make apparent that, contrarily to 
what one might infer from [EJ, it is not necessary to have flat space-time manifolds with (totally) compact spacelike 
sections M.3 for the topological induction of the reverberations and (or) the divergent behavior of the energy E(t) 
— the compactness of M3 in just one direction suffices. 

We begin the section ^ by emphasizing the importance of topological considerations in physical problems, and give 
an example that makes explicit how the overall electric charge of universes whose t = const spacelike sections are 
orientable and compact flat 3-manifolds is related to the topology. We also show that topological fragility can arise 
not only in the usual cosmological modelling, but also in the ordinary multiply-connected flat space-time manifolds 
we consider in the present paper. There we stress the origin of the topological reverberations we found, and show that 
it constitutes an example of topological fragility, which takes place without violation of any local physical law. The 
extent to which this paper generalizes ref. [EJ is discussed in details in section |y|. Finally the relation of our results 
with cosmology is also addressed. 



II. RADIATION REACTION EQUATION 

The flat space-time manifolds we shall be concerned with are decomposable into M4 = 1Z X M3 , where the real 
line 7Z represents a well defined global time of Minkowski space-time. Corresponding to the possible topologically 
distinct 3-manifolds M3 there exists a simply-connected covering manifold 1Z 3 , such that each manifold M3 is 
obtained from 1Z 3 by identifying points which are equivalent under the action of a discrete subgroup of isometries of 
the Euclidean space 1Z 3 without fixed points. In other words, each manifold M3 is obtained by forming the quotient 
space M3 = 1Z 3 /r, where T is a discrete group of isometries of 1Z 3 acting freely on 1Z 3 and referred to as the covering 
group of M.3 (for a complete classification of Euclidean 3-manifolds see Wolf 0[). 

The physical system we shall be concerned throughout this paper is represented by a pointlike harmonic oscillator 
coupled with a scalar field. This system has been treated by many people, since the work by Schwalb and Thirring Q 
in 1964, as a simplified model to study a few features of oscillating (radiating) phenomena || - |ll|. When the 
underlying manifold is 1Z 3 this system may serve as a model for an electric dipole coupled to electromagnetic radiation 
or an impurity atom interacting with acoustical waves [^) . It has also been used as simple model to parallel classical 
electrodynamics in the dipole approximation as discussed by Kampen Jl9| and Kramers p(| . It is worth mentioning 
that similar systems have also been used to mimic the basic properties, and thus to study the most relevant features of 
pulsating stellar systems J21J - |||] . The essential idea that permeates these latter papers is to work with simple models 
of oscillating fields interacting with oscillating sources in order to develop some analytical physical background to 
understand the relationship between gravitational waves and their sources (radiation emission and radiation reaction) . 

Consider that the harmonic oscillator of our physical system is located at an arbitrary point q £ M.3- Without loss 
of generality, one can assume that the covering map 7r : 1Z 3 — ► M3 maps the origin to q , namely 7r(0, 0, 0) = q. In 
what follows we shall denote by O q the orbit 7r _1 (g) of (0,0,0) under the action of T on 1Z 3 (see fig. 1). 

Now, the time evolution of the harmonic oscillator at q £ M3 can be obtained by studying the equivalent system 
on the covering manifold, which is formed by an infinite set of indistinguishable harmonic oscillators each one located 
at a point of the orbit O q , subject to identical sets of initial conditions and interacting with the scalar field ip in the 
same way as that of the original oscillator at q £ M3 (see fig. 1). Thus, the dynamics of our physical system in M.3 
can be derived from the following functional action for the equivalent system in the universal covering space 1Z 3 : 

S = S f + S + Si, (2.1) 

where Sf is a scalar field term, S corresponds to the oscillators, and Si is a coupling term between the scalar field 
and each harmonic oscillator. These three terms are, respectively, given by 

S f = \f d^x rT p(f , x) d v <p(t, x) , (2.2) 
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s ° = \ E / dt t Q} (*) - (*) 1 > ( 2 - 3 ) 

where i G (0, oo), <p(t,x) = <p(t,^x) is the massless scalar field, Qp{t) is the amplitude of the oscillator located at 
p G O q , overdot represents derivative with respect to the time t, u> is the angular frequency of each oscillator and 
A 7^ is the coupling constant. Hereafter, p is a non-negative normalized density function (J p — 1) with compact 
support (corresponding to the region of interaction) contained in the interior of a fundamental domain of A4$ and 
centered at the origin, and finally pp{x) = p(7 _1 x), where 7 G T and 7(0,0,0) = p [note that pg (x) = p(x)]. The 
step function 9 (t) is used to indicate that the interaction starts at t = 0. It should be noticed that the interaction 



term (2.4) clearly is not bounded from below, potentially permitting a limitless extraction of energy by the oscillator. 
However, as we shall discuss later in this article, it is the topology of the 3-space that "excites" this physically available 
mode of the system. 

A word of clarification is in order here. In dealing with the behavior of physical fields in the curved space-times 
of general relativity one has to consider that these fields are not only influenced by, but they also have gravitational 
effects, i.e., they change the geometry. In most ordinary situations, however, these back-reaction effects are small 
enough that can be neglected. In such cases the fields are treated as test fields in the corresponding curved space-time 
background. Had we considered our fields Q(t) and cp(t, x) as non-test fields in a non-static curved space-time we 
would have to take into account the back-reaction of the fields on the geometry, and so besides Si a new interaction 



term would have to be added to the action (2.1). However, our backgrounds are static flat space-times (endowed 
with any one of the 17 possible non-trivial flat topologies), and our system contains only test fields with a mutual 
interaction, thus there is no place for such a term in the action of our system. 



Varying the action (2.1) with respect to tp and Q one obtains the coupled equations of motion of the equivalent 
system, namely 

□ ip{t,x) = A PpW e ^ %w . ( 2 - 5 ) 

peo„ 



Qp (t) + uj z Qp (t) = A / d A x pp (x) 9 (t) ip(t, x) , (2.6) 



where □ denotes the d'Alembertian operator on 1Z 3 . Note that the sum in (2.5) is not an infinite sum. Actually 
owing to the disjoint compact supports of the functions pp (x) the right-hand side of that equation contains at most 
one summand which does not vanish. 



To obtain the radiation reaction equation from the above equations (2.5) and ( [2.6] ) we recall that the general 
solution of eq. (2.5), considered as an initial value problem, can be written in the form ip(t,x) = (pj(t,x) + <p H (t,x), 



where cp H (t,x) satisfies the corresponding homogeneous equation, and where ipj(t,x) indicates the solution of the 
inhomogeneous equation. From now on we shall assume, for simplicity, that the whole energy of the system is 



initially stored in the oscillator. Accordingly, we shall use as initial conditions for the equations (|2.5| ) and ( |2.6| ) that 
ip(0, x) = (p(0,x) — and that Q(Q) — a, Q(0) = /?, where a and f3 are real arbitrary constants. These conditions 
imply that f H {t, x) — 0, and therefore ip(t, x) — tpj(t, x). 

Now since the interaction begins at time t — 0, the Green function for the dAlembertian operator in ( |2.5| ) is clearly 
given byQ 



1 Clearly there are two equivalent approaches to build the Green function of the wave operator, which depend on whether it 
is considered on the quotient manifold M% = TZ 3 /T, or on the covering space 1Z 3 . On M3 , as it has been considered in 
it has to be invariant under the covering transformations. On 1Z 3 , however, which we are considering in this section, the 
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Hence, the solution for equation (2.5) can be formally written as 



f{t, x) 



A 

All 



reo„ 



pf(y) 



\x-y\ 



6{t-\x-y\) Q?{t-\x-y\) 



(2.7) 



(2.8) 



where we have used r as dummy index instead of p. 



Finally, inserting ( |2.8| ) into equation (2.6) one easily obtains the radiation reaction equations for the oscillators, 
namely 



? ? W+^ 2 %W = f £ / My ^ P0L e(t-\x-y\) QAt-\x-y\) 



F - y I 



(2.9) 



As the infinite set of identical harmonic oscillators (each one located at a point of O q ) are subject to the same set 
of physical constraints, the amplitudes Qp must evolve identically. Thus, the set of radiation reaction equations (|2.9| ) 
reduces to just one differential equation, namely 



em 



f -y I 



(2.10) 



which holds for each oscillator in O q C 7Z 3 , and whose solution Q(t) gives the time evolution for the amplitude of the 
harmonic oscillator in M3. 



One might think at first sight that the sums on the right hand side of eqs. (2.8) - (2.10) are divergent, and a 



regularisation scheme is needed to cope with the divergences. However, this is not true because for any finite time 
t > the sums in those equations are in fact finite owing to the disjoint compact supports of the functions pp (x) 
together with the cut-off effects of the step function 6(t). In other words, we do not have formal series to be carefully 
handled, but finite expressions. 

To deal with the pointlike coupling between the harmonic oscillator and the scalar field [ p (x) — > 5 3 (x) ] , we shall 
consider an infinite family of oscillator-field systems (n-system, for short), with the same coupling constant A =/= 0, 
and such that each element of the n-system is characterized by a bare angular frequency u> n and density function p n , 
satisfying the generalized Aichelburg-Beig condition 0], namely 



Air J \x-y\ 



n 2 > 



(2.11) 



where Q is the same constant for each member of the n-system, and {p n } is a 5-sequence in the sense of distribu- 
tions [f24j| . It should be noticed that the inequality ( 2.11 ) is a sufficient condition for having radiation damping of 
the energy E(t) in the pointlike coupling limit case when the oscillator-field system lies in the Minkowski space-time 
manifold M.4 = 1Z x 1Z 3 (see, for example, ref. M). 



Now, since eq. (2.1C) holds for each member of the n-system, using relation (2.11) one finds 

,(t)-0(T) Qn (r) 



Q„(i) + 2r / d i xd i yp n {x) Pn {y) 



\x-y\ 



n 2 Q n (t) = d 3 xd 3 y Pn{S) Pn ' A y ] 



feo„ 



\x-y\ 



(r) Q n (r) , 



(2.12) 



invariance under F must not be imposed. Note additionally that the covering transformations have been taken into account in 



eqs. (|2.S|) - (|2.6|) and in the subsequent equations (2.9) - ( 2.1S 
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where we have set r = t- \x-y\, 2Y = A 2 /4tt and O q = O q - {(0,0,0)}. 

In the limit n —> oo (see Appendix) we have that the following equation holds for all t > : 

oo 

Q{t) + 2TQ{t) + n 2 Q (t) = J2 °k 0(t - t k ) Q(t ~ t k ) , (2.13) 

fc=i 

where Ck = 2TDk/tk, Dk is the number of points r £ O q such that |r| = tk, the sequence {tk} is ordered by increasing 
values, and clearly t k — > oo when fc — ► oo. 



It should be noticed that the right-hand side of the radiation reaction equation (2.13) formally contains an infinite 



countable number of retarded terms of topological origin. However, due to the 9 function, at any given finite time t 



only a finite number of terms will contribute to the sum of ( 2.13 ). These terms are different for distinct multiply- 
connectedness of the 3-manifolds M3, and vanish for the simply-connected manifold 1Z 3 . Loosely speaking they 
give rise to "reflected" waves from the "topological walls" , and thus account for both the extrema of E(t) and 
the discontinuities of E(t) . In other words they are responsible for the topological reverberations. Clearly for the 
simply-connected manifold 7Z 3 , there is an exponential decay of E(t) and no reverberation takes place in this case. 

III. SOLUTION AND ASYMPTOTICAL BEHAVIOR 



In this section we shall present an exact recursive solution of the radiation reaction equation (2.13), which holds 



for any multiply-connected flat manifold To this end, we first consider the differential operator 

D = d 2 + 2Td + fl 2 , (3.1) 
which acts on real- valued piecewise smooth functions fit) according to 

Df = / + 2I7 + 2 / . (3.2) 
It is straightforward to show that if g(t) is some integrable function, then the solution of the equation 

Df = g (3.3) 

with initial conditions /(0) = /(0) = is 

fit) = er Tt f e TT h{t-r) g{r) dr , (3.4) 
Jo 

with /i 2 = Q 2 - T 2 and 

H~ Y sin fit if [i 2 > , 

h (t) = { t if /i 2 = , (3.5) 
v^ 1 sinhf t if /i 2 = — v 2 < . 



In what follows we shall discuss in detail how to find the solution of ( 2.15 ) for the underdamped case, in which fjr > 0. 
The solution for the remaining two cases, however, can be similarly worked out, and for the sake of brevity we will 
present later only the final results. 

Our purpose now is to find a real- valued function Q : [0, oo ) — + 1Z of class C 1 such that Q(0) — a and Q(0) = (3 



and satisfying (2.12), which in terms of the operator D can be rewritten as 

oo 

DQ(t) = J2 C k0(t-t k )Q(t-t k ) , (3.6) 



fe=i 



where the Ck are positive real numbers and {tk} is a strictly increasing sequence of positive real numbers (0 < t\ < 
ti < ■ ■ ■ < tk), such that tk — > oo when k — > oo. 



6 



Thus, for f(t) = Q{t) — So(t) and taking into account (3.3) and (3.4) the solution of (3.6) can be written as 



Q(t) = s (t) + -Y, C * *(* - **) ^ T(t ~ tk) 
M k=i 



e T sin fi(t — tk — t) Q(t) dr 



(3.7) 



where the function so must satisfy the homogeneous equation D so(t) = with the conditions so(0) =oe, so(0) = /?. 
Clearly one can rewrite the solution (p/fl) in the form 



(3.8) 



Q(t) = « (t) + £ c fe (* - **) <2i (* - **) , 



k=l 



where, from (3.3) and (3.4), the function 



1 r 

Q 1 (t) = - e _r * / e rr sin/i(i - r) Q(r) dr 
M Jo 



(3.9) 



fulfils the differential equation DQi (t) = Q(t) with the initial conditions Qi (0) = Qi(0) = 0. 

The solution (|3.8|) is given in an implicit form. An explicit form can be obtained, though. Indeed, as we are 
concerned only with real-valued functions defined for t £ [0, oo) , if we recursively define functions s„ by 



D s n (t) = s n -i(t) for n > 1 , 
Ds o (t) = 0, 

with the conditions s„(0) = s n (0) = for n > 1 , and functions Q n (t) by 

DQ n {t) = Q n - X {t) for »>1, 

Qo (t) = Q(t) , 

with Q n (0) = Q n (0) =0 for n > 1 , then we can show by induction that 

oo 

Qn («) = « B (t) + C * °( t - **) On+1 (* ~ *k) • 
fe=l 



(3.10) 
(3.11) 



(3.12) 
(3.13) 



(3.14) 



Notice that the first of these equations (for n = 0) is nothing but equation (3. 



An explicit form for the solution of (3J3) can now be obtained as follows. Equations ( 3.10 ) and ( 3.11 ) can be 
integrated to give 



Sq (t) = e ( A sin fit + B cos fit) , 



1 - f 

Sn(i) = — e rt / e TT sin fi(t — r) s„_i(r) dr for ri > 1 
A* Jo 



(3.15) 
(3.16) 



where 4 = (/3 + ar)//i,5 = «, and where obviously we have used (3J2) - (3ij). On the other hand, if we define 

go(*) = s (t) , 



'12 



(t) — 2_, Ck 2 o(t — t kl — tk 2 ) s 2 (t — t kl —tk 2 ) , 



fei ,fc 2 



(3.17) 



fci,fe 2> — ,ki \j=l / 3=1 3=1 
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then the explicit form of the solution of equation (p]q) is simply given by 



Q(t) = ]T fc(t) . 
=0 



(3.18) 



ake 



It should be noticed that since the above solution is ultimately given in terms of sq , equations ( 3.1 5| ) - (3.1 
apparent that Q(t) is completely determined by the pair of initial conditions [Q(0) = a , Q(0) — /3] as one would 
have expected from the outset. 



Another important point to be stressed regarding the above solution is that the right-hand side of ( |3.18|) is not an 
infinite sum as it appears at first sight. Actually, for a given finite time t the above sum will clearly be truncated 
by the presence of the step function [0(t) = for t < 0, 0(t) = 1 for t > 0], i.e. only those qi(t) with i <t j t\ will 



contribute to the sum (3.18j). Moreover, each qi(t) as given by (3.17) is itself a finite sum, again due to the cut-off 



effects of the step function 9{t) . 

Before proceeding to the discussion of the asymptotical behavior of Q(t) we mention that the solution of ( 2.1 3| ) for 
the critically damped (/i 2 = 0) and overdamped (/i 2 = — v 2 < 0) cases is again given by (|I§) with ( p7| ), but now 
with the functions Sq and s n given by 



*o(t) 

Sn(t) 



e~ Ft (A+Bt) for p 2 = , 
e~ rt (A sinh vt + B cosh vt) for p 2 



V < 



,-rt 



h(t — t) s„_i(t) dr for n > 1 , 



(3.19) 
(3.20) 

(3.21) 



where A and B are real constants, which are determined by the initial conditions, and the functions h(t) are given by 
equations (|3.5|). 



Now, for any non-trivial flat topology the asymptotic behavior of the solution Q(t) of equation (3.6) can be figured 
out by a procedure similar to that used in [pj. Indeed, if one ad hoc assumes an exponential asymptotical behavior 



for Q of the form Q(t) = 7 exp (at) with a and 7 real constants, then in the limit t — > 00 equation (3^) reduces to 



2 r a + n 2 



k=l 



k e 



(3.22) 



This equation has only one real solution for a. Indeed, the right-hand side of ( |3.22 ) is a positive monotone decreasing 
function r (a) with lim CT ^o T (a) = 00 and Imiu^oc r (a) = . Thus r (a) lies entirely in the first quadrant of the 
plane and crosses it from the top- left to the bottom-right (see fig. 2). Now, since T > then for a given pair (r, 51) 
the left-hand side of (3.22) is a parabola curved upwards with vertex at a = —T. Therefore, it always intersects the 
curve for r (a) in just one point, which is in the first quadrant (see fig. 2). This amounts to saying that there is only 



one real number a = b > solution to equation ( 3.22 ). Hence the amplitude Q(t), and consequently the energy of 
the oscillator E(t) = \ [ Q 2 (t) + £1 2 Q 2 {t)], both exhibit an asymptotical exponentially divergent behavior (confirmed, 
within the limits of accuracy of the plots, by figures 3 and 4 for the special case of an overdamped oscillator discussed 
below). One can sum up by stating that the compactification in just one direction is sufficient to topologically induce 
a asymptotic divergent behavior of E(t). 

A question which naturally arises here is whether the total energy of our system is conserved. As our system 
contains only mutually interacting test fields Q(t) and ip(t, x) in static flat background, the energy of the overall 
system must be conserved. Indeed, for t > we have 6(i) = 1 and the action ( pTl] ) - (2.4) does not depend explicitly 
on the time. Thus, for any well-behaved density function p(x) the total energy of our system £(t) clearly is conserved, 
which squares with the fact that our system is a test system in a static background. For the sake of completeness, in 
what follows we shall show in more details how the total energy of our system £ (t) is finite and conserved for all time 
t 6 (0, 00) even in the pointlike limit case p (x) —* S 3 (x) . We shall also qualitatively discuss how the balance of £(t) 
can be handled. 



From the functional action (2.1) - (2.4) the energy £(t) of our system can clearly be written in the form 
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8{t)=£ f {t) + £ {t) + £i{t) , 



(3.23) 



where 



d 3 x{<p 2 (t,x) + [V<p (t,x)} 2 }, 



£o(t) = -[Q 2 (t)+iu 2 Q 2 (*)] , 

£i(t) = -A Q(t) I d 3 x p(x)p(t,x) , 
Jv 



(3.24) 
(3.25) 
(3.26) 



for t £ (0, oo). Here V denotes a fundamental domain of As we have mentioned above for any time t > the 
energy £{€) is conserved since the corresponding Hamiltonian does not depend explicitly on the time. 

From equation ( |2.8| ) the interaction component £i(t), for which the field <p(t,x) is evaluated at x = 0, i.e. in the 
location of the pointlike oscillator, clearly presents a divergent behavior of the form |af| _1 in the neighborhood of x = 
(see below for details). Now, to deal with the balance of the total energy £(t), we introduce the divergent quantity 



. A y l 

A = — nm — 

47T r^0+ r 

in terms of which a finite energy Ei(t) can be defined for all < t < 00, namely 

£,(f)=B,(t)-AQ 2 (i). 



(3.27) 



(3.28) 



Thus, for each time t in that interval, clearly one can extract from £i(t) an infinite amount of negative energy to 
obtain a regularized finite negative energy Ei(t). 

one finds that in the neighborhood of x 



The energy £f(t) as given by ( 3.24 ) also diverges. In fact, from (2 
the scalar field ip(t, x) takes the form 



ip(t, x) 



47T Ixl 



+ (n.d.) 







(3.29) 



where (n.d.) denotes the non-diverging terms for any t S (0, 00). So, in the neighborhood of x = we have |V(/?| 



I (A/47r) Q(t) I |x| 2 |, and the integral corresponding to this term in eq. ( 3.24 ), namely \ J v 4Trr 2 dr IV^I 2 , has a 
divergent contribution of the form AQ 2 (t)/ 2 when r = \x\ — > 0. Now since the contribution of the term ip 2 to the 



integral on the right-hand side of (3.24) is finite for all t £ (0, 00), one can similarly introduce a regularized finite 
energy E f (t) by 



£ f (t)=E f (t) + -AQ 2 (t) 



(3.30) 



which is clearly finite for all finite time t > 0. 



The energy component corresponding to oscillator £ (t) as given by (3.25) can be dealt with by using the frequency 

S 3 (x) , namely using ui 2 



n given by eq. (2.11) in the limit case p (x 
frequency one can define 



£ (t) =E a (t) + - AQ 2 (t) , 



fi + A. Indeed, with this renormalized 



(3.31) 



where the regularized energy E Q (t) = | [ Q 2 (t) + fl 2 Q 2 (t)] is finite for all finite time t > 0. The energy of the oscillator 
E D (t) has also been denoted simply by E(t) in many places of this paper. 

Finally, we note that according to eqs. ( |3.28|) , ( |3.30| ) and ( 3.31 ), the divergent terms of the form AQ 2 (t) cancel out 



in the sum £j(i) +£ a (t). Thus, from ( 3.23 ) one obtains that the total energy of the system £(t) given by ( 3.25 ) 

reduces to 



£(t) = E f (t) + E (t) + Ei(t) , 



(3.32) 
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which is constant and finite for all time t > 0. 

The fact that the total energy of the system £ is finite and conserved for all time t > is undoubtedly the most 
important physical requisite our system obeys. However, it is the fact that the interaction energy is not bounded from 
below which makes possible and understandable the asymptotical divergent behavior of the energy of the oscillator 
and of the scalar field, since they both can "extract" energy from the interaction term. Remarkably this latent degree 
of freedom of our system is excited only when it lies in the multiply-connected manifolds, since when it is in 1Z 3 there 
is a damping of energy of the oscillator. In other words, the topology (multiply-connectedness) excites this available 
"physical mode" of our system. This remarkable feature of this system has not been perceived since the sixties, when 
it was first examined by Schwalb and Thirring (see, for example, refs. || - H]). 

It is certainly desirable that quantum physical systems have the interaction energy bounded from below. However, 
at a classical level this is not an imperative. When one considers, for example, a system of two isolated pointlikc 
charges with opposite signs in a flat space (vacuum medium), the smaller is the distance d between the charges the 
greater is the absolute value of potential energy of the system. Even if one assumes that the charges have already been 
renormalized, in the limit d — > the potential energy formally diverges. This sort of divergences are formal, since the 
limit d — >■ cannot be attained in practice. Furthermore, one ought to bear in mind that the classical theory used to 
infer this divergence does not hold for every small d (compared with, e.g., the Bohr radius). Actually, the distance 
d is assumed to be larger than the Bohr radius since otherwise the self-fields of the two interacting charges would 
be rather distorted, and one should expect much more involved physical effects. Similarly in our system the energy 
of the oscillator E Q (t) is finite for all finite t > 0, but formally diverges in the limit t — > oo. This limit case may 
perhaps be excluded on physical grounds. However, the most remarkable point here is that a well-behaved physical 
system when it lies in the ordinary simply-connected Euclidean manifold 7Z 3 , will exhibit an unexpected behaviour 
(reverberation pattern followed by a growth of E(t)) when it is in any possible flat manifold with non-trivial topology, 
without violating any local physical law. 



IV. CASE STUDY 



In this section we shall focus our attention on the time evolution of the harmonic oscillator in a flat space-time 
whose t = const spacelike section is the simplest multiply-connected Euclidean 3-manifold, namely A4^ = 1Z 2 x S 1 . 
In this case, denoting by a the distance between two equivalent points in the covering manifold TZ 3 , one can easily 
find that Ck = iT/tk, where tk — ka (k stands for a positive integer). Hence, from ( 3.17 ) and for i > one obtains 
the solution 



Qi(t) 



' k 1 ,k 2 ,...,k i 



iU 



9(t 



i=i 



Si (t — a kj 



(4.1) 



which together with (3.18) and the first eq. (3.17) give the behavior of the amplitude Q(t) and of the energy E{t) — 
\ [ Q 2 (t) + £! 2 Q 2 (t)] for the following different types of oscillator: (i) the underdamped ( n 2 > 0) for which so and s n 
are given by (3.15) and (3.16); (ii) the critically damped (a 2 =0) where (3.1£) and (3.21) furnish sq and s n ; (iii) the 



overdamped (u 2 < 0) for which s and s„ are given by (3.20) and (3.21). 

As far as the asymptotical behavior of Q(t) for this special case (M.3 = 7Z 2 x S 1 ) is concerned, the right-hand side 



series of ( 3.22 ) can easily be evaluated and equation ( 3.22 ) becomes 



2 T a + W 



4T 



ln(l 



(4.2) 



Clearly, for a specific oscillator, i.e. for a given pair (r, f2), this equation can be numerically solved for a by using, 
for example, a computer algebra system such as Maple |2^,|2^]. For an overdamped oscillator with r = 6 and £1 = 5, 
which we shall use in our numeric calculations, taking a — 1 and using Maple one easily obtains the approximate 
solution a = b = 0.35 . 



10 



In the remainder of this section we shall discuss the figures 3 and 4 and make some further remarks. We begin 
by emphasizing that we have chosen to discuss the dynamical behavior of both the amplitude and the energy of the 
oscillator in a space-time whose spacelike section is M% — 1Z 2 x S 1 because this is the simplest multiply-connected 
flat manifold (since it is compact in just one direction) having therefore the lowest degree of connectedness for a 
flat 3-manifold |5J]. This amounts to saying that the dynamical effects exhibited by the physical quantities in this 
manifold will also be present in any other multiply-connected flat 3-manifold. As a matter of fact, the greater is the 
degree of connectedness of the more reinforced will be these dynamical effets of topological nature. The limiting 
simply-connected case 1Z 3 will manifest no sign of such effects. A second point of general order which is worth noting 
is that we have decided to plot the graphs for an overdamped oscillator because we knew from the outset that the 
amplitude Q(t) and the energy E(t) for this oscillator, when in TZ 3 , exhibit an exponential decay with no relative 
extrema. This choice amounts to freezing the degrees of oscillations, and therefore rules out any sort of resonance, 
since clearly there is no frequency to be coupled with to resonate. Besides, our choice also makes easier and neater 
the comparison of the dynamical behavior of the oscillator in the simply and multiply connected cases. Taking into 
account these considerations we have plotted the figures 3 and 4 for an overdamped harmonic oscillator with r = 6 , 
51 = 5 and with initial conditions Q(0) = 1 and Q(0) = 0, in a flat space-time whose t — const spacelike section 



is endowed with the topology 1Z 2 x S 1 , and where we have taken the distance a = 1 in eqs. (4.1) and (4.2). The 
analysis of these figures is given in what follows. 

Figure 3 shows the graph for the amplitude Q(t) which exhibits relative maxima and minima followed by a growth 
of Q(t). This time evolution of the amplitude sharply contrasts with the behavior of the amplitude Q(t) for the same 
oscillator in TZ 3 , where neither relative extrema nor growth of Q(t) take place. This figure also contains the graph of 
the function exp(— bt) Q(t) , which within the accuracy of numeric calculations performed for the plots confirms an 
exponential form for the divergent behavior suggested by the graph of Q(t) and shown to indeed take place in the 



previous section. In this graph we have used b — 0.35, which is the root of eq. (4.2) for T = 6 , £1 = 5 and a = 1. 

Figure 4 shows the behavior of the energy E(t) which exhibits a few relative minima and maxima followed by an 
asymptotic divergent behavior. This divergent behavior is clearly in agreement with the behavior of Q(t) when t — > 00 
shown in figure 3. The time evolution of the energy for the oscillator in this multiply-connected manifold is in striking 
contrast with the behavior of the energy E(t) in the simply-connected case 7Z 3 , where there is an exponential damping 
of E(t) with no relative extrema. Figure 4 also contains the graph of the derivative E(t), exhibiting discontinuities 
at each t — tk — k a (k stands for a positive integer), which again do not occur in the simply-connected case 1Z 3 . The 
graph of E(t) also reveals that E(t) has a finite number of extrema. 

It should be also noticed that the derivative E = Q(Q + fl 2 Q) of the energy function exhibits discontinuities not 
only in the specific case handled in this section, but also in the general underlying topological setting of this work. 
Indeed, from equation ( 2.1 3| ) one easily obtains 



E(t) = 2TQ(t) 



00 „ 

J2~± e(t-tfc) Q(t-t k ) - Q(t) 

k=l k 



(4.3) 



From this equation one has that the discontinuities occur at t = t k , i.e. they come about each time a new term 
Q(t — tk) is taken into acount in the right-hand side of eq. ( 2.13| ). 



The relative extrema and asymptotic divergent behavior of E(t) as well as the discontinuities of E(t) , which occur 
only in the multiply-connected flat manifolds, are due to the retarded terms of topological origin in the evolution 
equation ( [2.13 ). These topological terms are different for distinct connectedness of the t = const section M.3 of the 



flat space-time manifold, and vanish for the simply-connected case 7Z 3 . They are ultimately responsible for those 
features (extrema, discontinuities and asymptotic behavior) of E(t). In brief, as it has been made apparent for the 
overdamped oscillator, both the suppression of the radiation damping and the reverberation pattern of E(t) in those 
multiply-connected flat space-times are of purely topological origin, and on the other hand, evince that our system is 
topologically fragile O] , 
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V. CONCLUDING REMARKS 



Since the physical laws are usually expressed in terms of local differential equations the topological considerations 
may not be prominent at first sight. Nevertheless they are often necessary in many problems of physics. In general 
relativity and cosmology, which handles the dynamic and the global structure of space-time manifolds, the study of 
topological features is even more significant and acquires a dynamical meaning in a sense. 

Topological considerations are necessary in many other situations in physics. When one considers, for example, the 
electric fields produced by bounded sources, one quite often chooses as boundary condition that the field vanishes at 
spatial infinity. This choice is possible and even convenient if the t = const spacelike section of the Minkowski space- 
time A^4 is the ordinary simply-connected Euclidean manifold 1Z 3 . However, as far as multiply-connected spacelike 
three-spaces M3 are concerned the choice of satisfactory boundary conditions for this problem is not as simple as that. 
So, for example, in the exam of the electric field produced by an isolated point charge in any one of the six possible 
orientable compact (multiply-connected) topologically distinct Euclidean 3-manifolds M3, clearly this condition at 
infinity cannot be imposed. Actually, it is easy to show by assuming the Gauss's law that one cannot have a net 
electric charge in these compact flat manifolds, making clear that the overall electric charge in these manifolds is 
related to the topology [0]. This example evinces that changes in the assumption of the space-time topology may 
induce important physical consequences even in the case of static flat space-time manifolds. This type of sensitivity, 
which is also present in cosmological modelling [jlTj , has been referred to as topological fragility and can clearly occur 
without violation of any local physical law jl 17|. As a matter of fact, in the above example the topological fragility 
arises when we impose the local validity of an ordinary law of classical electromagnetism to a system which lies in a 
flat multiply-connected manifold. 

In the same context, in this work we have studied the role played by multiply-connectedness in the time evolution of 
the energy E(t) of a radiating system that lies in static flat space-time manifolds M4 whose t = const spacelike sections 
A^3 are compact in at least one spatial direction. So, it may even have the lowest degree of multiply-connectedness. 
This topological setting is general enough to encompass the entire set of possible classes of multiply-connected flat 
3-manifolds discussed, for example, by Wolf (2). We have shown that the behavior of the radiating energy E(t) of the 
oscillator changes remarkably from exponential damping, when the underlying 3-manifold is 7£ 3 , to a reverberation 
behavior for E(t) when the spacelike t = const sections of Ad 4 are any possible multiply-connected flat 3-manifold. 
Thus, our study shows that the topological fragilities can arise not only in the usual cosmological modelling [fj"7| , 
but also in ordinary static flat space-time manifolds as long as they are multiply connected. As we have emphasized 
in section |l| and III this striking behavior of E(t) occurs with no violation of any physical law, and has a purely 
topological origin. 

Although the physical system we have investigated here is the same system discussed in 0, the present article 
generalizes the results of that paper in several respects. Firstly, the underlying topological setting is much more 
general since we did not assume that our spacelike 3-manifolds Ad^ are compact and orientable. Actually we only 
require that they are multiply-connected, which means that Ad^ can be endowed with any one of the 17 multiply- 
connected flat topologies ||], including, of course, the orientable compact 3-manifolds dealt with in Q|. Secondly, 
contrarily to the numerical integration performed in jl[] , here we have obtained a closed exact solution of the evolution 
equation for the harmonic oscillator in the above-mentioned general topological setting. Thirdly, in the case study we 
discussed in details in section IV we have found the (exact) time evolution of our system in a specific non-compact 
(although multiply-connected) manifold, while in Q the time behavior of our system was (numerically) treated only 
for six compact orientable background manifolds. Finally in the present article we have also discussed the conservation 
and the balance of the total energy of our physical system, making apparent that in both papers the total energy is 
finite and conserved for any time t > 0. 

Before closing it is worth mentioning that to the extent that we have studied the dynamical behavior of a test 
system in static flat FRW space-time backgrounds our results have a moderate, though suggestive and well-founded, 
relation to cosmology. Besides, it should be noticed that the net role played by topology is better singled out in this 
static case, where the dynamical degrees of freedom have been frozen. On the other hand, it appears that, apart from 
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the inflationary expansion and perhaps a few other cases, the restriction we have made will not be decisive for the 
pattern of the behavior radiating energy E(t) in a number of flat expanding locally homogeneous and isotropic FRW 
cosmological models. Finally, it is worth noting that since geometry constrains, but does not dictate the topology of 
space-time manifolds, in cosmological modelling one often is confronted with the basic question of what topologies are 
physically acceptable for a given space-time manifold. An approach to this problem is to study the possible physical 
(observational) consequences of adopting particular topologies for the space-time (for a fair list of references on this 
approach see the review article by Lachieze-Rey and Luminet fig}] ). This paper constitutes an example of this relevant 
approach. 
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APPENDIX: POINTLIKE COUPLING LIMIT CASE 



Our aim in this appendix is to show how one can derive from equation (2.12), which holds for an arbitrary density 



function p n (x) of a ^-sequence, the radiation reaction equation ( 2. IS ) for the limiting case of pointlike coupling 
between the scalar field and the oscillator [ p [x) — > S 3 (x) ] . 
Let us first rewrite the equation ( [2.12 ) in the following form 



Q n (t) + 2T I n (t) + n 2 Q n (t) = 2T J2 Jn,f(t), (A.l) 



r e i 



where we have set 



TM frf ^ Qn(t)-e(r)Q n (r) 

In(t)= Id x d y p n (x) p n (y) p — , (A.2) 

J \x-y\ 

J n ^(t) = [ d 3 x d 3 y Pn{ f} Pn f® 9(r) Q n (r) , (A.3) 

J \x-y\ 

with T = t-\x*-y\,2T = A 2 /4tt and O q = O q - {(0,0,0)}. 

As we mentioned earlier the pointlike limit case can be achieved by taking the limit n — > oo of the above n-system 
equation. To this end we firstly remind that as we are dealing with density functions p n (x) whose supports are 
connected sets, the following reworded version of the mean value theorem can be stated: let f(x,y) be an arbitrary 
continuous function defined in 2 = supp (p„) x supp (p n .r) , then there exists a point (xq , yo) G X such that 

f(xo, Vo) = J d 3 x d 3 y p n (x) p n ^ (y) f(x, y) (A. 4) 

holds. 

Let now R n be the radius of supp (p n ) and consider a point f € O qi then by the mean value theorem there exists 
a sequence of pair of points (xo. n , yo.n) G T- such that for each t > we have the sequence of numbers 

if t<\r\-2R n , 

Jn,F (*) = { unknown if | r\ - 2 R n < t < \ r\ + 2 R n , (A. 5) 

Qn(t-e n )/e n if \f\+2Rn<t, 
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where e n = \ XQ, n — yo. n | . Since R n — > when n — > oo , the mid subinterval shrinks to the point t = \f\, then it is 
irrelevant for our purpose the knowledge of J n p(t) in that region. Furthermore, since e n — » | r | we have 



lim J n j (t) = 



8(t-\r\)Q(t-\r\) 



(A.6) 



The limit of I n (t) can be obtained as follows. For a given fixed value t > 0, there exists a natural number no 
such that for all n > no we have t > 2R n . Thus by the mean value theorem there also exists a pair of points 
xo,n , Vo, n € supp (p n ) such that 



In(t) = 



Qn(t) - Qn(t - e n ) 



(A.7) 



with e„ = | 2?o, n — yb,n | — * when n — > oo . Since Q n {t) has continuous first derivative we can expand it about the 
fixed value t up to first order to obtain 



Therefore the following limit: 



I n (t) = Q n (t) + 0(e n ) 



lim = I n (t) = Q{t) 



(A. 



(A.9) 



holds for alH > 0. 

Finally, using equations (AT), ( A.6 ) and ( A.S) one obtains the radiation reaction equation ( 2.13 ) for the pointlike 
limit case p n (%) ^ S 3 (x) . 
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FIG. 1. Two-dimensional schematic figure of a pointlike harmonic oscillator at a point q £ Mi and its equivalent system on 
the covering manifold 1Z 3 , which is formed by an infinite set of identical oscillators each one located at a point of the discrete 
orbit 7T _1 (q) — O q C IZ 3 . All oscillators interact identically with the scalar field y> and are subject to the same set of initial 
conditions. 

FIG. 2. Schematic representation of the curves corresponding to the left- and right-hand side of equation (3.22), whose 
intersection in just one point evinces that this equation has a single real positive solution a — b. 

FIG. 3. Graph of the amplitude Q(t) for the overdamped harmonic oscillator with T = 6 and Q — 5 in a flat space-time 
whose t — const spacelike section M3 is endowed with the topology 7Z 2 x S 1 , and where we have taken a = 1 in eq. (4.1), and 
as initial conditions Q(0) = 1 and Q(0) = 0. There are a few relative maxima and minima followed by a growth of Q(t). This 
time evolution sharply contrasts with the behavior of the amplitude Q(t) for this overdamped oscillator when M3 = 1Z 3 , where 
there are neither relative extrema nor growth of Q(t). This figure also contains the graph of the function exp(— bt) Q(t) , which 
within the accuracy of the plot confirms an exponential form for the divergent behavior suggested by the graph of Q(t). To 
have this last graph we have used b = 0.35, which is the approximate root of eq. (4.2). 
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FIG. 4. Behavior of the energy E(t) for the overdamped harmonic oscillator with T = 6 and Q. — 5 in a flat space-time 
whose t — const spacelike sections M3 have the topology 1Z 2 x S 1 , and where we have taken a = 1 in eq. (4.1), and as initial 
conditions Q(0) = 1 and Q(0) = 0. The graph for E(t) exhibits a few relative minima and maxima followed by an asymptotic 
divergent behavior. This divergent conduct is clearly conform to the behavior of Q(t) when t —* 00 . The time evolution of 
the energy for the overdamped harmonic oscillator in this multiply-connected manifold is in striking contrast with the behavior 
of the energy E(t) in the simply-connected case 1Z 3 , where there is an exponential damping of E(t) with no relative extrema. 
This figure also contains the graph of the derivative E(t), exhibiting finite discontinuities at each t = tk = k a (k stands for a 
positive integer), which again do not occur in the simply-connected case 1Z 3 . The graph of E(t) also reveals that E(t) has a 
finite number of extrema. The relative minima and maxima of the energy E(t), the discontinuities of the derivative E(t), and 
the asymptotic divergent behavior, which take place in multiply-connected flat manifolds, all have a topological origin. 
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